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ABSTRACT 

Tne PERT/COST system, now in current use, requires seweral simplify 
ing assumptions in order to permit the use of analytic wethods for deter- 
mining predictions of critical path activities and project duration. As 
a result of these assumptions, the predictions may involve significant 
error, and will always be optimistic. The Monte Carlo technique of net- 
work calculation does not require these assumptions and, hence, is 
capable of yielding more accurate predictions and providing more useful 
information. This technique is discussed in detail. A more flexible 
probability density model for activity times is introduced, and a ree 
source allocation technique based on the probability that an activity 
will be on the critical path is developed, Finally, an application of 


PERT network theory to military operational pianaing is described. 
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NOTATION 


Meaning 
Indicates a random variable 


Time required to complete an activity not yet 
completed 


A specific value of the random variable t* 
Expected value, or mean, of t* 


Time allowed by schedule for completion of an 
activity 


Elapsed time till the occurrence of an event, or 
until completion of an activity 


Expected value, or mean, of T* 


Elapsed time at which an event is scheduled to 
occur 


Latest elapsed time an event can occur and not 
delay the scheduled completion of the project, or 
the latest time an activity can be completed and 
not delay the project 


Expected value, or mean, of T,* 


Earliest elapsed time an activity/event can be 
Scheduled for completion 


Latest elapsed time at which an activity/event can 
be scheduled without causing a change in t, for 
some subsequent activity 

Slack for an activity/event, = T¡* - T* 


Expected value, or mean, cf K* 


Standard deviation of the random varíable in 
parentheses 


Variance of the random varíable ín parentheses 


Duration of the entire project 
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NOTATION (cont'd) 


Meaning 


Elapsed time to the Directed Completion Date for 
the project 


Probability an activity/event will be on the 
critical path. The Criticality Index. 


Probability am activity/event will be completed 
on schedule 


Probability that the project will be completed 
on or before the DCD 


Directed Completion Date for the entire project 
Optimistic estimate for t*, the smallest amount 
of time the activity can reasonably be expected 


to consume 


Most likely estimate for t*. The elapsed time at 
which the activity is most likely to occur. 


Pessimistic estimate for t*. The greatest amount 
of time the activity can be expected to consume. 


The standard Beta variable 


The mode of x* 
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m Introduction. 

The PERT system was devised by a study group under the direction of 
the Special Projects Office, Bureau of Naval Weapons, Navy Department, 
in order to develop a methodology for providing the management of the 
Fleet Ballistic Missile (Polaris) program with an information reduction 
system for program monitoring and evaluation. By this system, management 
was to be continuously apprised of progress to date for the program as a 
whole, and be furnished valid predictions as to outlook toward accom- 
plishing program objectives. [2] 

The system, as developed, was a giant stride in management tech- 
nology, and is credited with a major contribution toward the rather 
phenomenal success of the FBM development program. Not only did the 
PERT system accomplish the objectives already mentioned, but it allowed 
management to predict those activities in the development praject whose 
completion time would have a direct effect on the overall project dura- 
tion. These critical activities could then receive appropriate mana- 
gerial attention. 

The PERT system has since been employed extensively in industry, 
and has been made a standard procedure for monitoring research and 
development projects under the cognizance of the Federal government. 
Industry acceptance of the system has been, in general, enthusiastic and 
widespread. 

The system has been recently extended to include a cost control 
feature, PERT/COST, which has now been designated the standardized sys- 
tem for monitoring federally sponsored R&D projects. 


Another system, the Critical Path Method, has paralleled PERT in 





development. It is similar in many respects to PERT, lacking PERT’s 
stochastic representation of time, but incorporating a correlation be- 
tween time and cost, enabling management to schedule optimally. The 
Department of Defense and NASA have incorporated the desirable features 
of both systems into PERT/COST. 

The success of PERT and Pert-type systems has stimulated the inter- 
est of professionals in the fields of mathematics, management science, 
and operations research. The professional journals have featured many 
papers on the relative merits of the systems, proposals for refinement 
and extension of the systems, consolidation and integration of the sys- 
tems, and the like. It is the purpose of this paper to consolidate many 
of these ideas, and to propose some refinements and extensions to the 
System, designed to extend its usefulness and increase the accuracy of 
predictions about the projects, and to generate specific recommendations 
for control of time and cost. 

We shall first briefly describe the basic PERT system, critically 
discuss certain features of the system, and tender proposals by which 
the basic system may be improved. In particular, we shall discuss the 
Monte Carlo method of calculation of the PERT network, and show how this 
method overcomes many of the deficiencies inherent in the analytic method. 

Next we shall describe the standard PERT/COST system currently being 
implemented by DOD and NASA. We shall follow this with our proposal for 
an integrated project control, scheduling, and resource allocation tech- 


nique. 
Finally, we shall discuss an extension of the PERT network approach 
to the problem of planning projects involving major future decisions, 


such as military campaigns. 





2. Conclusions and Recommendations 

As a result of this study, it is concluded that; 

(a) The assumptions employed in the currently standard analytic 
PERT network calculations may induce significant error in the results. 

(b) The standard analytic technique for network calculation fails 
to provide a sufficiently definitive measure for activity criticality. 

(c) A more flexible model for activity completion time, providing a 
range of variance to allow for varying uncertainty in time estimates, is 
needed, 

(d) The Monte Carlo technique does not depend on the assumptions 
which induce the inaccuracy in the analytic technique, hence is capable 
of yielding results limited in accuracy only by the validity of the 
estimates for activity time. In addition, this method yields the proba- 
bility that a given activity will be on the critical path, and other 
useful information not obtained by the analytic technique. The Monte 
Carlo technique can easily provide the flexibility discussed in (c), 
above. 

(e) The Monte Carlo technique requires increased computer running 
time, but the increase is not prohibitive, and is well justified by the 
improvement in quality of the results obtained. 

(f) A resource allocation technique, based on the probability that 
an actívity is on the critical path, is feasible, and provides a reason- 
able method of optimum resource allocation. 

(g) The network representation of the complex interrelationships 
between activities in a project is a technique which may well be applied 


to the task of planning military operations. 





It is recommended that the Monte Carlo technique for PERT network 
calculations be employed in lieu of the analytic technique. A method of 


utilizing this technique is discussed in detail in Appendix I. 





3. The Basic PERT System 

The basic PERT system is a means of defining the relationship be- 
tween the various activities comprising a project, and for estimating 
the time required to complete each activity and the project as a whole. 
Use of the system enables realistic schedules to be imposed on the 
project, and assists in controlling the execution of the project. 

By a project we shall mean a collection of tasks, or activities, 
each a necessary step toward the achievement of some final objective. 
Subsets, called paths, of these activities, are dependent sequentially, 
that is, any activity in the path, other than the beginning activity, 
may not be commenced until its predecessor has been completed. Each 
activity in the path, except the beginning activity, has a predecessor 
in the path. An example of a project in this sense may be as simple as 
the construction of a house, where the activities are tasks such as lay- 
ing foundation, erecting subfloor, installing rough plumbing, erecting 
framing, installing rough wiring, etc. Or a project may be as complex 
as development of the Fleet Ballistic Missile system, with many hundreds 
of activities involved. 

A network of directed arcs is a very convenient way of representing 
a project. The activities are represented by directed lime segments, 
called arcs, terminating at nodes, called events. (See Figure 1.) The 
time required to complete an activity is associated with the length of 
the arc representing it. Events represent points in time. Several 
activities may terminate at one event. If so, the event is said to occur 
when the last of these terminating activities is completed. At that 


point in time, other activities, whose commencement is contingent upon 
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Figure 1 


A Network Representation of a Project 


the event, may commence. These activities are represented by arcs di- 
rected outward from the event. In the house construction example, the 
installation of the interior wallboard may commence only after the rough 


plumbing and rough wiring are completed. 
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Wall board 


Figure 2 


Two Independent Activities Incident to an Event 


A project can always be represented by a network with one beginning 


event and one end event. Fictitious, or dummy, activities may be drawn 





in the network in order to achieve this format. Such dummy activities 
may consume zero time, or may represent a predetermined amount of sched- 
uled "dead" time. It is common practice for the network to be drawn 
with beginning event on the left, and with the end event, representing 
the achievement of the final objective, at the extreme right of the net- 
work. 

A path may be defined as an unbroken chain of events, with origin 
at the beginning event and terminus at the end event. The network con- 
sists of many such paths, overlapping, paralleling, and crossing each 
other. 

Activities generally consume time, The exact amount of time, t, 
which will be consumed in completing any one activity is mot generally 
known in advance. If one looked at the history of a completed project, 
one could discover the exact amount of time, t, required to complete any 
given activity, and could label each event with its actual time of occur- 
rence. The longest path, with respect to time, from the beginning event 
to the end event could then be found. This path is called the critical 
path. The activities on this path are called critical activities. The 
length of the critical path, the sum, Dt, of the completion times for 
all activities on the path, is the total elapsed time from the beginning 
to the end of the project, and is called the duration, .نا‎ It is apparent 
that a small change in t for one of the activities on the critical path 
causes the same amount of change in D, but a small change in t for an 
activity not on the critical path would cause no an in D. Hence it 
is important that project management be able to predict the critical 


path in order to predict D, and ın order to allocate resources, and fix 





schedules, most effectively. 

Any path, other than the critical path, is shorter than the critical 
path. The difference in lengths of the critical path and amy other path 
is called the slack, K, for the non-critical path. An activity or event 
may be Located on several paths simultaneously. The smallest slack 
associated with any path which includes a given activity or event is the 
value of slack associated with that actívíty/event. The activities/ 
events on the critical path have zero slack. Any single activity may be 
delayed an amount less thàn or equal to its slack without affecting the 
duration of the project. 

The term, slack, may have a slightly different meaning if some au- 
thority has imposed a Directed Completion Date, DCD, on the end event of 


the project. If so, we define the Scheduled Duration, D., as the elapsed 


3? 
time from the beginning event to the DCD. Then slack for the critical 
path is (D, - D). Slack for any other path is (D, - path length). With 
this definition, slack for any activity may be either positive, negative, 
or zero, but all activities on the critical path will have the same 
slack, which will be less than that for any other activity im the net- 
work. 

The PERT system is a method by which predictions of the critical 
path, and slack for any activity, may be made with sufficient accuracy 
to enable management to operate more effectively. 

The Stochastic Network Model. In the planning stage, before à proj- 
ect has commenced, the actual time, t, which will be required to complete 


any given activity, is an unknown quantity. It is unlikely that this 


time could be predicted exactly. We shall call such uncertain quantities 





random variables, and symbolize them with an asterisk, as t, The corre- 
sponding unstarred symbol will represent a specific value of the variabie 
It is possible to set limits, within which t* is almost certain to fall. 
It is also possible to estimate the shape of the probability density 
function governing the random variable, t*. The probability density 
function for t* is a curve whose ordinate in a region is a measure of 

the likelihood that t* will fall in that region. To be more precise, 

the area under the curve in a region is the probability that t* will 
occur in that region. in Figure 3 the ratio of the shaded area to the 
total area under the curve is the probability that t* will fall between 


t and t5. 


1 
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Figure 3 


A Probability Density Function 


The actual probability density function for t* is difficult, if not 
impossible to obtain, and will be unique for any given activity performed 
in a given time period. Theoretically, if the activity were performed 


repeatedly, under the same conditions, with no learning taking place, 





then a histogram showing the number of completions occurring in a time, 
t t At, plotted against t in discrete segments, would approximate the 
probability density function. Since this procedure is not possible, we 
must ke satisfied with the best estimate we can formulate, using avail- 
able knowledge of the nature of the activity. 

In order to estimate the probability density function of t* for a 
given activity, the most qualified supervisor in charge of the activity 
makes three estimates of the completion time for the activity. One 
estimate, Ó, an optimistic estimate, is the smallest amount of time the 
activity may be reasonably expected to consume. Another, M, is the 
amount cf time the activity is most likely to require. Finally, he 
makes a pessimistic estimate, P, which is the greatest amount of time 
the activity can be expected to consume, barring completely unforeseen 
circumstances, P and @ form the upper and lower bounds, respectively, 
of the variation of t*, and M is the mode of the demsity function. The 
probability density function is assumed to have the value zero at P and 
6, and to reach a single peak at M, 

Any of several well known mathematical functions exhibit the prop- 
erties we have just assumed for the distribution of t*. The Beta 
function, 

f(t) = k(t-ay* (b-t) (1) 
with parameters a,b,«, and Y , was chosen as the model for the PERT sys- 
tem activity time because it could be specifíed to fit the constraints 
of the estimates. The three estimates 9, M, and P, specified three of 
the parameters of the Beta function. The remaining parameter was speci- 


fied by arbitrarily requiring the standard deviation, O (t*), to equal 


10 





22 , One-sixth of the range is a frequently used estimator for the 
standard deviation of unimodal frequency distributions, hence was con- 
sidered to be a reasonable assumption. For example, the standard normal 
distribution, truncated att 2.66, has its standard deviation equal to 

1/6 the range. [6] As a result of the above requirements, one of the 
parameters, say A , becomes a function of the relative position of M in 
the range between @ and P. Determination of this parameter requires the 
solution of a cubic equation, (see Appendix [). The exact determination 
of the mean, or expected value of t*, symbolized by t, , is a tedious 
process, since te is a function of a,b, A and Y. The following reasonably 


close linear approximation for the mean was adopted as a standard: 


E [t*] = te s LME (2) 


The probability density function thus assumed has a shape similar 
to that shown in Figure 4. This distribution will be referred to subse- 


quently as the PERT Beta distribution. 


qo 





Figure 4 


A Typical PERT Beta Density Function 


A more detailed discussion of the PERT Beta distribution is included in 


Appendix I. 
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For purposes of predicting the critical path, the three estimates 
are combined by means of the weighted average formula, (2), for the mean. 
The expected value, or mean, ta, is used as a deterministic value of t* 


in the prediction of the critical path, t, lies between M and the mid- 


e 
point of ($,P), and is separated from M by 1/3 the distance between M 
and the midpoint. In this manner, account is taken of the usual tendency 
for estimators to be optimistic about placing M. The fundamental reason, 
however, for using tẹ instead of M as an estimator for t* lies in proba- 
bility theory, rather than an effort to correct for this supposed bias, 
In order to determine Slack, two new random variables, T* and 7 
are defined. T* is associated with a given event, and is the elapsed 
time between the beginning event and the occurrence of the given event. 
T* is the sum of the random variables, t*, for the activities on the 
longest path to the event. At this point in the development of the 


a 


event will be the path whose mean length is greatest. The mean length of 


SSS = SaaS Se ee Se ee SS سس‎ 


a path is the sum of the means, t of the activities on the path. Then, 


e 
by the above assumption, the mean, Tg, of T*, is the largest of the meau 
lengths of all paths to the event. Operating on the basis of the same 
assumption, the further assumption is made that T* is approximately nor- 
mally distributed. This assumption is justified since we assume we know 
the longest path to the event, and since this path is the sum of n random 
variables, the Central Limit Theorem of probability theory asserts its 
tendency toward normality as n increases. Due to the assumed independ- 


ence of the activities along this path, the variance, G ب‎ of T* is 


the sum of the variances, deter), of the activities on the path to the 


12 





event. The distribution of T* is now completely specified by Tp and 
S CTH) ín view of the normality of T*. 

We may associate a random variable, T*, with an activity as well as 
an event. We define T* for an activity as the elapsed time until the 
activity is completed. Note that T* for an activity is not necessarily 
equal to T* for the event following the activity. T* for that event is 
egual to the maximum of the T's for the activities incident to it, that 
is, terminating at the event. 

The random variable, T,*, for an event is defined as the latest 
time the event can occur and not delay the scheduled completion of the 
project. TS is calculated by subtracting from Dg the sum of the t*'s 
for all the activities on the longest path from the event to the end 
event. Again the assumption is made that this Longest path will be the 
path whose mean length, DUM is longest. Then the mean value of T,%, 
written T,, = (Dg ~ 2٢ where the sum is the maximum of all such sums 
on paths between the given event and the end event. The variance of 
a o *(T)*), = c ^(t*), along the path with greatest mean length. 

A 1andom varíable, T,*, may also be associated with an activity. 
T,* for an activity is defined as the latest elapsed time since the be- 
ginning event at which the activity must be completed in order not to 
delay the scheduled completion of the project. Notice that T,* for any 
activity is the same as the B2 for the event following the activity. 

It is apparent that the sum of the two paths used to determine T* 
and Tj,* for a given event/activity is the longest path from the begín- 
ning to the end which contains that event/activity. From our definition 


of slack, we obtain the following relationship for slack, K*; 
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Ke Dg ˆ [m و‎ 1,+)] : ۳ (3j 
Then the mean value, Kg» Of slack is the difference, T, Te Notice 
that this expression for Kp depends upon the assumptions made ín deter- 
mining T; and Tg: 

Ihe Probability of Meeting a Scheduled Completion Date. We have 
developed the random variable, T*, whose distribution we have specified, 
by means of simplifying assumptions, to be normally distributed, with 
mean Tp and variance, O ery. The probability, Po, that an event will 
occur before some scheduled completion date is simply ۳ ]1* > 11 , where 
Tg is the elapsed time from the beginning event to the scheduled comple- 
tion date. [In particular, the probability, Pp. that the project will be 
completed before the directed date ís p [D* < D]. Where the duration, 

D*, is a random variable equal to T* for the end event. These probabili- 
ties can be obtained from a table of values for the normal random 
variable. 

The Predicted Critical Path is that path whose mean value of slack, 
Kp, is smallest. All events or activities on this path will have the 
same value for Kg, which will be less than the Kg for all other events / 
activities. 

This completes the description of the basic PERT system which is 
currently in general use. Tn practice, the network representing the 
project is drawn, events and activities are assigned identification num- 
bers, and are assigned to specific organizations for responsibility. 

The supervisors responsible for completion of each activity then submit 
the three estimates, ó, M, and P. Several computer programs have been 


written which take the activity number, preceding and succeeding event 


14 





numbers, and the three time estimates as inputs, then perform the net- 
work calculations. [5,7,15,20] The output lists the activities in any 
predetermined ordering, giving Teo Ty» Kp» the probability, Pp: of meet- 
ing a DCD, dry, ta, and g"^ ces). The report may list the predicted 


critical path, and several other paths in increasing order of K Man- 


ES 
agement may take appropriate action by reallocation of resources if the 
probability of meeting the overall project DCD is umacceptably small. 
In order to "buy time" additional resources are allocated to, and addi- 
tional managerial attention devoted to those activities on the predicted 
critical path, and to those whose mean slack is relatively small.  Re- 
sources may be shifted from paths with greater Kg to paths with lesser 
Kp. In some cases reorganization of the network is required, placing 
more activities in parallel rather than in series paths, or eliminating 
activities whose products are not entirely essential to the project. 

The results of original PERT computations may be the basis for 
establishing an overall project schedule. This schedule specifies sched- 
uled completion times, t,, for activities, and scheduled times of occur- 


rence, Ta, for certain events. 


g? 
While the project is in progress, revised inputs are submitted 

periodically, showing percent of completion, and revised estimates. In 

the case of a partially completed activity, these estimates should be 

considerably more reliable than those submitted before work was commenced. 

These inputs are fed to the computer, and another network calculation is 

run, as before, except that completed activities are noted and the actual 

completion time is used in the network rather than an estimated time, 

The results of these runs may necessitate schedule changes, and/or any 

of the actions previously mentioned. 


L5 





2 Critique of the Basic PERT System 

The system which we have just described tends to identify and call 
attention to a single path, the Predicted Critical Path. This path, 
though actually uncertain, is determimistically established by virtue of 
several assumptions, Im most networks, the activities on this "critical 
path" will have probabilities of actually being critical much less than 
unity. Furthermore, other activities, not on this path, may have proba- 
bilities of being critical larger than some activities om the "critical 
path". [18] Hence, the system, by focussing attention on the single path, 
may deemphasize activities which are actually more important. 

The assumptions which were made in determining Tp, and Ty lead to 
optimistic results in these figures, and in Kp. By reason of the same 
assumptions, the probability that an event is completed before a DCD is 
always optimistic. 

T* is the random variable representing the elapsed time to occur- 
rence of an event. Suppose there are n different paths to the event. 
Let the random variables, p1*, p>*, “aes be the lengths of these n 
paths. Now T* is the maximum of these n variables, but the maximum of a 
set of random variables is not normally distributed, even if the indi- 
vidual variables themselves are normal. In particular, T* does not have 
the same distríbution as the p* vith the largest mean, nor does T* have 
the same mean as this p*. Furthermore, the p,* are neither independent 
nor identically distríbuted, hence the calculation of the distribution 
of the maximum of this set is so complicated as to render it infeasible 
by analytíc methods. 


In general, the probability that am event occurs before some 
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specific scheduled date is the probability that none of the possible 
paths overrun the scheduled date. By the PERT assumptions, we consid- 
ered only one path, namely that with the greatest mean Length. Ob- 
viously, if there are several parallel paths with nearly the same mean 
length, then by considering only the probability that one of the paths 
does not overrun, we are being quite optimistic in our result. 

For example, suppose there are three parallel and independent paths, 
each approximately equally likely to be the critical path, and each hav- 
ing a probability of being less than D, of about 0.50. Then Py the 
probability that the project is completed prior to the DCD, is the prob- 
ability that none of the path lengths exceed D4. Hence, Pp = (.5)° = 20% 
The PERT procedure would consider only the path with the largest mean 
length, and would calculate Pp = 0.50. For networks with several paral- 
lel paths of comparable length, the error would be quite large, as in 
the example. In networks with only one predominantly long path, the 
error would be negligible. 

In order to obtain a more general prediction of criticality than 
the single path prediction, we need to calculate the probability that 
each activity will be on the critical path. We shall call this proba- 


bility the Criticality Index, written, P Computation of Pc involves 


C 
the concept of determining the probability that a particular Bi will be 
the maximum of a set of n random variables, fp rh. 4S before, these وم‎ 
are not independent nor identically distributed, hence this computation 
is not feasible analytically. 16 is not, however, difficult to deter- 


mine Po by the Monte Carlo method. This method also readily yields 


Tg, Ty, and Pp without resorting to the assumptions which were so 
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troublesome, but necessary, in the analytic approach. <A detailed dis- 
cussion of this method, and suggestions for specific use of the criti- 
cality index will form a substantial portion of this study. 

The PERT Beta probability distribution was chosen as the distribu- 
tion for t* for an activity. Recall that im specifying this distribution 
the variance of the distribution was arbitrarily chosen to be 8 
By this assumption, we have specified the degree of uncertainty with 
which the estimator makes his three estimates. Yn other words, we are 
denying the difference between the case where the estimator may predict, 
with a high degree of confidence, that t will occur within a few time 
units of M, and the case where t has a higher Likelihood of occurring 
near the extremes of the range of t*. The probability distribution for 
the former case should have a shape similar to that shown in Figure 5, 
while the latter case would be better represented by Figure 6. It is 
not unreasonable to assume that, under certain conditions supervisors 
will have a basis for estimates with a narrow uncertainty band about the 
mode, even though the range of possible values of t* might be quite 
large. Under other conditions, the supervisor might be unable to place 
M with an appreciable degree of confidence. in the interest of more 
accuracy in predicting the system behavior, perhaps a choice of density 
functions should be available to the estimator to help describe his con- 
fidence in his estimates. This would not require a knowledge of proba- 
bility theory on the part of the supervisor. He could simply be asked 
to rate his uncertainty as to the relative position of M on a numerical 
scale, say with three descriptive choices. His choice would lead to a 


corresponding probability distribution. Appendices Y, Ii, III and IV 
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describe means of providing this flexibility in choice of density func- 


tions with the desired degree of variance. 
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Figure 5 


A Density Function with Small Variance 





Figure 6 


A Density Function with Large Variance 
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p The Monte Carlo Technique 

We have seen that the completion time, t*, for any activity ín a 
proposed project is a random variable. As the project is performed, 
each activity is completed in some actual time, t. Only after comple- 
tion of the entire project can actual values of t, slack, and elapsed 
time, T, for events and activities be fixed, and the critical path de- 
termined, 

The Monte Carlo method calculates average values of these param- 
eters by repeatedly simulating the performance of the entire network, 
and statistically averaging the results. The simulation of the project 
is accomplished as follows: 

A random value of t for each activity in the network is drawn from 
che probability distribution for t* for the activity. These random 
values of t, called realizations of t*, are drawn in such à way that in 
repeated drawings from a specific probability distribution, the number 
of realizations falling within a segment of fixed length on the t axis 
is proportional to the height of the probability density curve over that 
segment. Each activity has its own unique probability distribution, de- 
termined by the three estimates, hence, for one simulation of the project 
a realization of t* drawn for an activity is a completion time which 
could have resulted from a normal performance of the activity. One simu- 
lation of the project consists of drawing a value of t for each activity 
of the network. Calculations of the network are then performed, using 
these values of t instead of the t, used in the analytical method. T and 
K for each activity are calculated, and then the critical path for that 


Y 


one simulation is determined. For each activity a tally is incremented 





if the activity was actually on the critical path Another tally ıs in- 
cremented for each scheduled event if the evemt was accomplished prior 
to its scheduled completion date. 

The above simulation is repeated many times and the values of t, T, 
and K for each activity are added to cumulative totals after each simula- 
tion of the project. After an adequate number of replications, these 
cumulative totals and tallies are divided by the number of replications, 
yrelding the average values, te, Te and Kp for each activity and the 
probability that the activity will be on the crítical path. For ali 
scheduled events, the probability that the event was completed prior to 
the scheduled date is also computed. The observed variance of these 
Statistics is also available if desired, as well as an approximate proba- 
bility distribution for project duration. 

The method of sampling the probability distributions of the activity 
times for realizations of t* is discussed in detail in Appendices I, Ll, 
III, and LV. Using the above techniques, it is not difficult to provide 
a choice of probability distributions for t* in order to account for 
varying degrees of uncertainty im estimates, as was discussed earlier. 

The computer codes designed to implement the basic PERT system may 
be modified, without great difficulty, to perform the Monte Carlo calcu- 
lation, The author recently modified the North American Aviation PERT 
Mark III code in this manner. This program, designed for the IBM 7094, 
will accommodate a network of 1500 actívities without requiríng external 
memory units. 3 

The Monte Carlo technique is not subject to the pitfalls of inde- 


pendence, assumptions of normalcy due to the Central Limit Theorem, and 
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assumption of a particular Longest path, etc., which pisgue the analytic 
computation. The inaccuracies inherent in the Monte Carlo technique are 
«ue to non-randomness in the random number generator, and perhaps, faíl- 
wre to perform a sufficient number of replications to permit the law of 
large numbers to obtain. Excellent random number generators are avail- 
able, providing results well within the Limit of accuracy imposed by the 
time estimates. Statistical methods exist by wnich to determine the 
number of replications required in order to obtain the desired degree of 
confidence in the results. (17 | 

The system is Limited in accuracy only by the human error in estimat- 
ing activity times, and in approximating the correct probability دم‎ 
bution for t*. 

The Monte Carlo technique requires more computer running time than 
does the analytic method. An efficiently coded routine requires approxi - 
mately 750 microseconds to draw a random number and perform the calcula- 
tions necessary to transform ít into a realization of t*. For details 
of the method employed, see Appendix I. The approximate increase in 
PERT Program time required, if the Monte Carlo technique is used in lieu 
of the standard analytic method, may be estimated by the product: 

(No. of activitıes)(No. of replications)(750 x 107€ sec). 

For example, a network consisting of 1000 activities would require 
12.5 minutes of additional computer running time if a Monte Carlo calcu- 
lation of 1000 replications were used instead of the standard technique. 
The advantages inherent in this technique would seem to render tais 1n- 
crease in running time negiigible. 


a 
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6. The Extension to PERT/COST ¡10 

An extension of the basic PERT system, PERT/COST, was adopted 1 July 
1962 as the basic standard for management time and cost control systems 
by DOD and NASA. PERT/COST extends the original PERT system to provide 
cost estimation and cost monitoring features in addition to the procedure 
for prediction and control of time performance. Two additional optional 
features, a Time-Cost supplement, and a Resource Allocation supplement 
were also promulgated. 

PERT/COST adds to the basic system the following capabilities: 

2. Initial estimation of project cost, broken down by work pack- 

ages, sub-systems, etc. 

b. Provision of a consolidated estimate of the requirements for 
manpower and other resources as am aid to scheduling and pro- 
curement. 

C Provision of in progress cost reports, showing current cost 
status and revised cost estimates in comparison with budgeted 
expenditures and contract estimates. 

The optional Time-Cost supplement provides a method of estimating 
three feasible schedules for project completion, with their estimated 
costs and an evaluation of the technical risks associated with each 
schedule. These schedules are designed to provide completion in accord- 
ance with three criteria: (1) the most efficient use of time and re- 
sources; (2) completion by the DCD; (3) earliest feasible completion. 

The Resource Allocation supplement outlines a procedure for sched- 
uling a project to meet a DCD in an optimal manner with respect to costs, 


and utilization of manpower. 
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The concept of operation of PERT/COST is essentially as follows 

A. Project breakdown. The project is broken down into activities 
and work packages, and graphically represented by a network. For the 
purposes of cost reporting, estimating, and accounting, several activi- 
ties may be combined into a "work package", when the activities involved 
are small, and detailed cost reporting and estimating of these activities 
would involve unnecessary expense. In gemeral, a work package does not 
comprise a closed network. That is, it may not be integrated and repre- 
sented by a single activity due to the presence of events interior to the 
work package which are connected by activities to other work packages. 
In this case, the work package must be broken down into its component 
activities and events for the network time calculations. The cost esti- 
mating and expenditure reporting section of the System may utilize the 
entire work package as the smallest organizational division. 

Certain work packages may not appear on the project network. In 
general, these are the functions of management not directed toward the 
accomplishment of specific activities and events. For example, the 
accounting, purchasing, or management travel costs for the project would 
be difficult to represent as activities in the network; however, the 
cost of these work packages may be estimated and reported, and treated 
in the cost analysis portion of the PERT/COST system. 

B. Time calculations. Time estimates are made for each activity, 
and the network calculations performed. The project duration is com- 
pared with the DCD, and if necessary, any of the following actions taken 
in order to reduce the overall project duration: 


(a) Reorganize the network by increasing the parailelism of activi- 


ties. 
24 





(b) Commit additional resources to selected activities, or reallo- 

cate resources within the network. 

(c) Change or delete activities. 

C. Manpower loading. [9] A breakdown by skills of manpower required 
for an activity/work package may be included with the time and cost esti- 
mates. Computer programs now in existence assemble this information and 
prepare a report and display showing the project manpower requirements, 
by skills, plotted against time. 7ط‎ This information provides project 
management with a forecast of manpower needs, and a basis for reschedul- 
ing activities in order to best utilize available manpower, Scheduling 
a large projeet without considering the overall use of manpower skills 
leads to very uneconomical manpower use. One time period may require a 
manpower level for a certain skill far in excess of that available, while 
the foilowing time period makes little use of that skill. Im performing 
۰ necessary rescheduling, the time constraints of the network and 
criticalness of activities must be considered. Those activities least 
likely to be critical are the obvious candidates for rescheduling for 
the purpose of leveling manpower requirements. Computer programs now in 
existence, supplementing standard PERT routines, perform this reschedul- 
ing under the constraints of maximum availability of manpower and with 
consideration of slack. (7, 11) 

D. Scheduling the network. A scheduled duration, tg, 15 promul- 
gated for each activity. These times may be less than, equal to, or 
greater than te, depending on the a and management policy. Using 
the ts for each activity, the network is computed and values of Sp and 


Si, corresponding to Tp and و1‎ determined. SE represents the earliest 
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date on which an activity/event can be scheduled for completion. Sr 
represents the latest date on which an activity/event may be scheduied 
for conpletion without causing a schedule slippage for the proiect com- 
pletion. SE and S, are computed in exactly the same manner as were TE 
and و1‎ when t. is used as activity time instead of te. As a final step 
in scheduling, an elapsed time, Tg, is imposed 0و"‎ event. To is 
chosen to lie between SE and St» and may be translated into a firm sched- 
uled cempletion date for activities immediately preceding the event, and 
a beginning date for activities immediately following the event. 

E. Budgeting the project. After the schedule is finalized, cost 
estimates are prepared for each work package, and a budget prepared for 
allocating funds. 

F. Monitoring progress of the project. As the project progresses, 
reports are submitted from activities/work packages, indicating progress, 
new time estimates, funds committed to date, and new estimates for cost. 
This information forms inputs to the PERT/COST computer program, which 
performs the network calculations, and prepares reports showing new 
values for Tp, T, Kp, Pp» and the new critical path, as well as current 
expenditures in comparison with estimates and budget. The new projected 
cost curve is also generated and compared with estimates and budget. 
Armed with this information, management may reschedule, reallocate, re- 
budget, or take whatever action as may be appropriate to the situation. 

The Time-Cost Option Supplement. The purpose of this supplement is 
to standardize procedures by which a project manager may prepare three 
alternative project schedule proposals, and evaluate the technical risks 


associated with each. By technical risk, we mean the gamble of 
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performance, time, and cost incurred by departing from the best develop- 
ment technique in order to meet the schedule. 

First a plan is formulated designed to meet the project DCD, 5 
imposed by the controlling agency. 

Next, a Most Efficient plan is formulated, in which the best devel- 
opment techniques are used in order to reduce technical risks, and most 
efficient use is made of manpower and resources, resulting in greater 
duration for the project. Total cost is generally lower for this plan. 

Finally, a Shortest Time plan is proposed. The purpose of such a 
plan is to explore the feasibility of scheduling the project for a dura- 
tion shorter than that proposed by the contracting agency in order to 
benefit from any strategic benefits which might accrue from such early 
completion. In order to arrive at the shortest time plan, the schedule 
is compressed by allocation of additional resources, paralleling activi- 
ties which should normally go in sequence, eliminating activities, and 
changing technical approaches to the problems of design and construction. 
All such actions serve to increase the technical risks involved, and most 
of them increase cost. 

The three options are presented to the contracting agency with esti- 
mates of cost and time, and with a comparison and evaluation of the tech- 
nical risks associated with each plan. 

The PERT/COST system is the primary tool for obtaining the time and 
cost estimates needed to prepare the three schedule proposals. 

The Resource Allocation Supplement. This supplement provides a pro- 
cedure by which project managers may schedule the project in the optimal 


manner with regard to the time-cost trade-off. The concept of operation 
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is as follows. 

For any activity there may exist several feasible schedules, involv- 
ing different time-cost relationships. (See Figure 7.) Efficient use 
cf manpower and existing resources and machinery may result in schedule 
A. Further time extension to B might result in higher costs due to the 
effect of fixed costs. By use of overtime, hiring additional manpower, 


or purchase of additional machinery or space, time-cost combination C or 


D might be achieved, 


— — —e Ms 
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time: 


Figure 7 


Time-Cost Plans for an Activity 


Schedule A would be the most economical for this activity, hence 
would be the time estimated for the activity for the first network cal- 
culation. In a similar manner, most economical estimates are obtained, 
with corresponding cost estimates, from the other activities. The net- 


work calculation is then performed, and the project duration, D, compared 
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with Dg. Ef D>Dg, selected activities on the critical path must be 
Shortened by resorting to a higher cost-shorter time work schedule 

For each activity on the critical path, determine the increase in 
cost, divided by the decrease in time that would occur in moving to each 
shorter time-cost point. This is the slope of the time cost curve for 
the corresponding time reduction, The activity with the lowest value of 
slope is then chosen as the activity to rescheduie. 

If this time reduction in the critical path exceeds the difference 
in slack between the critical path and the smallest slack not on the 
critical path, it becomes necessary to recompute the network and deter- 
mine a new critical path. The above process is repeated as often as 
necessary until DZ De. 

Slack paths are then reexamined to see if any activities may be ex- 
tended to lower cost points without going critical. If fixed costs were 
not a part of the activity estimates, notice may be taken of the fact tha: 
the most economical point to operate is usually to the left of the mıni- 
mum cost point with respect only to direct costs. After determining the 
most economical plan with respect to direct costs, a test is made to de- 
termine if further time reduction below the DCD would result in lower 
overall costs. This is done by adding on the fixed costs for the project 
and attempting further reduction in time. If further time reduction re- 
sults ín lower overall cost, then further time reduction is performed 
until total costs reach a minimum. 

The schedule is then adjusted as necessary to reach an optimum lev- 
eling of manpower and other resources within the existing constraints 


and considering total cost. 
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7, A Method for Resource Allocation, using the Criticality Index 

In Section 4, we pointed out the errors in the concept of a deter- 
ministic prediction of the critical path. We also showed the optimistic 
bias inherent in the assumption that the longest path to an event will 
be the path whose mean length is greatest. The Monte Carlo technique 
does not depend on this assumption, hence, with this technique we can 


calculate Pc, Tp, and P_ with accuracy limited only by the supervisor's 


D 
estimates and our approximation of the true probability distributions of 


t* for the activities. We shall now develop a system for optimum resource 


liccation utilizing the criticality index obtained from the Monte Carlo 


a 


calculation, 

{n order to use this method, we must require that a work package be 
a closed network. A closed network is a network with a single beginning 
event and a single end event. Activities exterior to the closed network 
may not be incident to, or emanate from, events within the closed network 
other than the beginning or end event. (See Figure 8.) In this example, 
if any activities not included in the work package should be incident to, 
or emanate from events B, C, or D, then the network would not be closed. 
This same restriction is implicit in the resource allocation supplement 
of the standard PERT/COST system, although not stated in the PERT/COST 
document. Notice that the closed network may be integrated and repre- 
sented by a single activity connecting events A and E. Some PERT com- 
puter programs now in existence are capable of integrating such sub- 
networks for simplification of the network for the benefit of higher 


9 9 


levels of management. /15! 
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Figure 8 


A Closed Network 


We next require that the supervisor in charge of each work package 
submit as many time-cost plans as are feasible for his work package com- 
pletion. These plans are represented by points on the time-cost curve ín 
the same manner as previously díscussed in connection with the Resource 
allocation supplement. No points to the right of the minimum cost point 
are submitted. The time-cost information is represented as a set of n 
pairs, (time,cost), in descending order of time and corresponding ascend- 
ing order of cost. For example: 

14 weeks, $40,000 
12 weeks, 60,000 
10 weeks, 90,000 

The time corresponding to the minimum cost point is taken as the 
most Likely estimate, M, of activity completion time. The additional 
time estimates, @ and P, are also submitted for the minimum cost plan. 
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These three estimates and the corresponding cost are the basis of the 
preliminary calculation of the network. From this calculation, the preb- 
ability, Pp, of meeting the project DCD, is determined. Recall that this 
probability is considerably more valid as a result of the Monte Cario 
technique. 

If Pp is unacceptably small, the network must be compressed by re- 
designing the network layout and/or allocating more resources to selected 
work packages in order to achíeve shorter completion time. If, after 
whatever network redesign is considered feasible, the network must be 
further compressed, the resource allocation phase of the program is per- 
formed. This phase consists of the following logical steps: 

l. For each time-cost plan of each activity having positive Pc, 


the following figure of merit is computed. 


ETB = Pc کے‎ (4) 


where At is the reduction in time achieved by using the T-C plan being 
considered, over that used in the previous network calculation. AC is 
the corresponding increase in cost. ETB may be considered as the expected 
amount of time which will be "bought" per unit cost by selecting the T-C 
plan being considered, 

2. Form 3-tuples consisting of each ETB, its corresponding At, 
and its activity identification. Order this list in descending order of 
ETB. 

3. Select from the top of the list, a sufficient number of reduced 
time plans so that the sum of these A t's does not exceed a predetermined 
percentage, p, of D. The best value for p may be determined experimen- 
tally. The selection is made from the top of the list, but only one 
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plan, that with the largest At, is retained for any one activity. 
4. For the activities chosen for time reduction, the time corre- 
sponding to the new T-C plan is taken as the new value of M. The esti- 


mates for P and Ø are revised downward by multiplying by the ratio, 


M (new) 
M (old) 


would be to require P and Ø estimates for each T-C plan.) 


(This procedure may be subject to question; another method 


5. With the new time and cost estimates for the accelerated activi- 
ties, the network is again calculated, and Pp determined. All activities 


will have a new P The mean duration, and Pp are again evaluated. if 


c 
still unacceptable, the process may be continued until either an accept- 
able Pp is attained, or the budgetary limit is reached, The schedule is 
finalized on the basis of the tíme plans used in the last network calcu- 
lation, after consideration of manpower and other resource leveling. 

By the above method, resources are allocated where, probabilistic- 
ally, they can be expected to contribute most effectively to shortening 
project duration. The requirement for several feasible time-cost plans 
may serve the additional purpose of forcing line supervisory personnel 
to consider the time-cost relationship more carefully, allowing operation 
at the most efficient point on the time-cost curve whenever feasible. 

An optimum value for p can be determined with a few experimental 
runs of the program. p should be chosen as large as possible in order 
to decrease the number of network calculations required, but small enough 
so that successive calculations produce a smooth compression of the ner- 
work. Operation of the program will tend to decrease the differences 
between the ETB's for the activities, which is equivalent to operating 


the activities at the same level of expected marginal utility. This is 
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a well known principle of optimization. 

The inputs to the system, three time estimates and a T-C function 
specified by several 2-tuples, are not complicated, and are compatibie 
with the concept of progress and cost monitoring envisioned in the PERT/ 
COST system. 

A flow diagram of the logic used in this resource allocation method 


is shown in Appendix V. 
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5 A Military Application of the PERT Network Approach 

One of the most outstanding features of the PERT approach is the 
graphical portrayal of the project by means of the network. By this de- 
vice, relationships between the various tasks may be clearly visualized. 
The insight thus gained by the managers of the project, into these com- 
plex relationships, enables them to function with greatly increased 
effectiveness in their managerial capacity. 

A military operation is in many ways similar to a development proj- 
ect. The operation may be divided into a number of tasks, with many 
interrelationships existing between individual tasks. Hence, the net- 
work may serve the same useful function in this application as in indus- 
trial development. 

The phases of military planning known as "The Development of the 
Plan," and "Supervision of the Planned Action," are most readily bene- 
fited by the network approach. In the formalized military planning 
sequence, the Development of the Plan phase has been preceded by the 
"Estimate of the Situation," in which the assigned mission has been 
studied, possible outcomes based on enemy capabilities and alternative 
own courses of action have been analyzed in a game theoretic matrix, 
and a decision has been reached regarding a general course of action to 
be pursued. 

The problem now involves developing a complete, detailed plan for 
the operation, including organizing the available forces and assigning to 
each task unit the appropriate tasks which make up the general plan. At 
this stage, detailed planning is done in which training, acquisition of 


intelligence, movement, communications, logistics, anc battle action are 
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ali considered. 

The plan begins with an initíal concept, which may be represented 
by a relatively simple network, whose tasks are stated in broad terms of 
accomplishment. As the planning becomes more detailed, these initia! 
tasks may be subdivided into networks of lower-order tasks. This pha 
may be done at lower echelons of command, upon receipt by them of thei 
superior's directive. 

The network will aid in establishing a firm schedule for the oper 


tion, in much the same manner as is done in industrial projects. t may 


also be a valuable aid for conducting briefings for subordinates, since 


it will enable them to better appreciate the relationship of their own 
tasks to the whole operation. 

During the operation, the network may serve as a display, on which 
current progress is indicated as reported. The effects of delays and 
failures may be more readily evaluated and corrected. 

The initial concept of a plan to execute an amphibious landing in 
enemy-held territory might be represented by the network shown in Fis 

A network may also be used to lay out the basic strategy of a pro- 
tracted military campaign. In thís usage, tasks would be represented by 
subsidiary operations, such as: "Seize Island B," or "Cut supply line 
XY.' The procedure is similar to that employed in the operation just dis- 
cussed, except that it takes place over an extended period of time. in 
many cases, proposed actions are tentative, or the selection of alterná- 
tive actions may depend on future developments, such as the success 
failure of a previous endeavor, or on the subsequent enemy reaction je 


may indicate this need for decision at a particular point in the network 
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by showing a "decision box'" in lieu of an event at that point. The deci- 
sion box has more than one task emanating from it; the decision, to be 
made at some later date, determining which task(s) to be performed. The 
network then proceeds to show the planned action to be taken for each 


alternative. Figure 10 is an example of a portion of such a strategy 


network. [13,14] 
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APPENDIX I 
Monte Carlo Techniques with the Beta Distribution 
The model for activity completion time, t*, chosen by the originators 
of PERT, is the Beta probability distribution, characterized by its den- 


sity function: [2] 


A K(t-0)" (P-t)Î , 6۵ < م ع‎ (5) 


= 0, elsewhere. 
We may reduce this to the standard form of the Beta distribution by the 


following transformation: 


x* = 5 (6) 





The probability density function of x* is: 
OK رک ۶ 0 سر‎ (7) 
The variable, x*, is equivalent to t*, where @ = O and P = l. 
The two parameters, @ and P, specify the upper and lower limits of 
the variation of t*. It is now desired to specify the remaining two 


parameters, & and or so that the mode of the density function occurs at 


M, and so that the distribution has the desired variance. 





Let m = Ls denote the mode of x*. By finding the root of £ 220) 


we obtain: 


A 
wo. ۶ (8) 


The variance of x* is given by: 


je: با وا‎ a 
(oL K-- 232 (oc * + 3) 





Now suppose we desire that the standard deviation, d (t*), be expressed 


oy (232, where d is an arbitrary constant. Then کل‎ (x*) . i. By 
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. 1 
requiring that (یں)2 ک‎ mn و‎ We obtain the following relationship in A 
and m. : 
+7 2 d'm*u Gap 0 
3 2 3 a 2 . / [| سب م‎ )m oc ` 
E 4 (d me —d^m E 7m)oc 7 ( (10) 


Now œ can be determined as a function of m and d by the above equation. 


With « determined, / is specified by solving equation (8): 


1 == (11) 
The mean of x* is given by: 
كن‎ F 
EN ST Xm 12 


Since X must be determined from the cubic equation, (10), for each value 
of m, the computation of E(x] is cumbersome. In the PERT Beta model, d 
was chosen as six. For this value of d, E[x] is approximately linear ín 


m, and may be approximated by the following relation: 





E [x] mm! (13) 


The transformation to t* yields: 


+ 4M P 
E [c] 4 £ AM FF : 


(14) 
Random Sampling from the Beta Distribution. There are several good 
random number generating routines which yield random numbers from the uni- 
form (0,1) distribution. We shall use the technique of the Probabiiity 
Integral Transformation to transform the random number, u, drawn from the 
uniform (0,1) distribution, to a corresponding sample from whatever other 


probability distribution we desire. ]3[ 


This transformation depends on the following theorem: 
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For any random variable, v*, having the probability distribution 


random variable, 


function, F,4, and the density function, f xo define the 
u* z Fs? that SE 


V 


|J £ Ct) dt (15)‏ کہ َو 


— 00 


mee EO EEE A ee ee EEE E ë ey e Ėė Á ےہ س‎ o 


Hence, we may use the inverse function, v = A to transform the 
sample, u, from the uniform (0,1) distribution to the corresponding sample, 
v, from the desired distribution. 


The distribution function: 


x 
3 1 | 
Fkx(3) = Ky | y (1-y) dy (16) 
0 

for the Beta variable, x*, must be calculated by numerical integration. 
In order to utilize the procedure given above for sampling from this dis- 
tribution, we may store the tabled distribution functions in the memory 
of the computer and use a table look-up routine to enter the table with 
the random sample, u, and determine the corresponding value for x. A 
complete distribution function must be stored for each of the increments 
of m, and for a limited number of values of d. 

The input card for an activity specifies the choice of d which best 
represents the estimator's uncertainty. Also specified are the three 
estimates, Ø, M, and P. With d and M specified, the proper table may be 
selected. u is produced by the random number generator. The correspond- 
ing value of x is obtained from the table, then transformed into a reali- 
zation of t* by the following transformation: 
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AA‏ لاجم 


ا ل .€ HAS‏ 








un = Da > -‏ هه سس 


t a (P Ona (175 
Tables I through V are tables of distribution functions for ti« 
standard Beta distribution, with d taking on values of 4 through 8, re: 
spectively. A complete function is given for each of 11 values of = 
AS an example of the use of the procedure given above, suppose (, M, 
and P were 12, 15, and 20 weeks, respectively, and the selection, d = 6 


-1 
is made. Then m = = : = .375. Suppose we draw the random number, 





.795, from the uniform random number generator. Entering Table [if for 
m= .3, with the argument, .795, we obtain x = .50. From the table for 
m= .4, Linear interpolation yields x = .572. Interpolating between 
these two values of x, we obtain x = .554 for m = .375. The correspond- 
ing realization of t* is then 8x +12 = 16.4, 

An efficiently coded routine designed to perform the above tabie 
searches and calculations, requires approximately 750 microseconds. 

Figures ii, 12, and 13 are families of curves of the probability 


density functions for x* corresponding to d = 4, 6, and 8, respectively. 
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M is the position of the mode 
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M is the position of the mode 
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Figure il 
Beta Density Functions for d = 4 
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Beta Density Functions for d = 6 
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Figure 13 
Beta Density Functions for d = 8 
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APPENDIX II 
A Simpler Method for Operating With the Beta Distribution 
Since the Monte Carlo method does not require that the variance of 
t* be calculated and applied in order to evaluate Pp, we may consider dis- 
pensing with the requirement that the variance of t* be constant over the 
range of M. This requirement is not necessarily a natural one. The vari- 
ance of the triangular distribution, for example, defined naturally by O, 


M, and P, is given by; 


2. 
سے‎ 2(t*) = 2 2-6 Z (P-M (M-Ó (18) 
6 18 
UB see that the variance for this distribution varies from a maximum of 
Em? 3(P-0 wi 
لس‎ with m at the extremes of its range, to 3 (Bo with m at the 


midpoint. The variance of the SinCos function, (Appendix IV), another 
function naturally specified by the three estimates, behaves in the same 
general manner. 

We may simplify the determination of the Beta function parameters, 
avoiding the necessity to solve the cubic equation, (10), by the follow- 
ing procedure: 


“ix the variance of the symmetric Beta function, (m = .5), by 


letting Ä : 
4 *(t* z 2201 (19) 


۱ d 7 
iM E و‎ k 


where d is an arbitrary positive constant = 12 0 


Tr Evo (20) 
d2 
m <=.) 
For the case where m = .5, equations (8), (9), and (20) yield the 


following: 
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2 
o E. 5 = d^ - 12 (21) 


2 
Define C = Lu (223 


Then, for any value of m in its range, let & and Y be determined as 


follows: 


(23) 


The remaining parameter may be determined by equation (8). The variance 
ef x*, as determined by equation (9), then varies with m as snown in 


Figure ll. 





Figure 4 


Change of Variance with Mode Position 


The mean, E (x* , may be calculated directly by equation (12). The 


following close linear approximations for E [x*], for several values of d, 
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rA 


show fue relative weights given to Ø, M, and P in the determination of 
expected value, hence give an intuitive feeling for the relationship be- 


Cween d and the amount of uncertainty with which m is located. 


A El (24) 


4 20 * M ۳ 
2 

6 Ó * 4M + P 
6 

8 Ó X 6M E 
8 


By this relatively simple procedure, the parameters of the Beta 
function for any desired values of d may be determined for the purpose of 
the Monte Carlo calculations. The Probability distribution functions 


needed for the transformation from u* to x*, may be extracted from a table 


e 


۶ the Incomplete Beta Function. [1] 
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APPENDIX III 


Use of the Triangular Probability Density in PERT 
Monte Carlo Calculations 


the triangular probability density function is naturally defined by 
the three time estimates, Ø, M, and P. Since the Distribution function of 
this density may be calculated by a simple formula, the transformation 
from a sample, u, drawn from a uniform (0,1) distribution, to a sample 
from the destred triangular density, may be accomplished by a single, 
simple calculation. When this procedure can be used, a considerable sav- 
ing in computer running time is effected. 

the technique is illustrated below for the random variable, y*, hav- 
ing the general triangular distribution with lower límit, 0, upper limit, 


P, atid mode, M. 





Figure 15 


Triangular Density Function 


The density function of y* is: 


Ó —y cM‏ نس 


2( 
(P-0) (M-0) 





$ ماله 7 
) > 


۱۱ 


Bere) > 2 
(P-B)(P-M) ۰ er 


- 0, elsewhere 
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/y 
fhe distribution function, Pye) = E yu (2) dz, is: 


Ww 00 


۱ E ۰ 
yO) = TOM) دج‎ 


IN 


M 


۱ )26( 
_ MeO, M- 2PM+ 2Py - y? , Me y &P 
۵۔۶‎ (P-0) (P -M) 


Now let u = Ey (Y) > and. solve for the inverse function, y = Fy (u), 


and we get: 


(27) 6ھ و ے 0 , Ø + Vu(P-0)(M-0)‏ 


< 
u 


P - vy(P-M)(P-Ø)(l-u) , Es £u ے‎ 1 


XE we draw a random number, u, from the uniform (0,1) distribution, 
and determine y by means of the function, (27), then y* has the desired 
triangular distribution. 

The mean, E(y*] , is given exactly by: 


g[y*] = d (28) 


The variance, g ?(y*), is given by: 
(yx) = ان‎ E کا‎ Cueca (29) 
The term on the right varies from 0, when M is at either extreme, to 


و 
seu , when M is at the midpoint.‏ 





APPENDIX IV 
The SiuCos Function 


Another Probability Function Suitable for PERT 
Monte Carlo Calculations 


Another probability function which is completely specified by its 
extremes and its mode is the following, which we shall call the SinCos 
function, for obvious reasons. 


The probability density function is: 


f a(z) "m" Tr (2-0) 


E ۰ سے‎ ٠ 
= zT ° 6۵ عء‎ 2 5 ۸ (30) 
= A TL(z-M Me = 7 

2)۳-0( ۰ 20۳-1 ' 
= و0‎ elsewhere 


Tre function has the following shape, where M can take any position be- 


tween $ and P. 


7 ہے 
* 
E‏ 
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Figure 16 


The SinCos Density Function 
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7 


the distribution function, F,,(2) = emer, is; 
La 
F,4(2) = 0, 2 > ۵ )31( 
. M-Ó h.cos Tp (2-9) 0 ج ے‎ ۷ 
۳-۵ 1 2(M-0) | ” 
< ۱۳ PM sin (2-M MZE P 
P-0 ۲ 9 2(P-M) ' 7 
= 1 z= P 


Setting u = F,4(2), and solving for z, we obtain: 


N 
8 


- 6 + 2(M-8) cos”! |. ے 0 , اسف‎ M-Ó 
TY (M-0) | P- 


OIS 


(32) 


X 
IN 


2(P-M) 0. -ı |P-8/ _ MOA ۲-۵ 
M + Tr sin EX 5-6 > 52 u £ 1 


The above function transforms a random sample, u, from the uniform (9,1) 
cistribution to the corresponding sample, z, from the SinCos distribution. 
The mean, El: of the SinCos distribution is: 


E2 = MDOE ODE (33) 


Considered as a weighted average, we see that the extremes, @ and P, re- 
ceive more weight than M. 


The variance, (٭ج)2 ک‎ is: 


T.‏ رف ود جم تی کی O‏ مت 


In a form for numerical comparison with the PERT Beta variance: 


2 
rer) s 2.06(2-8) - „364 (E0069) (35) 





We see that 6 *(2%) is 1.88 + .18K times the PERT Beta variance, P-6%, 
6 


where -1£ K< 1, depending on the relative position of M. 
Although, with the SinCos function, the transformation from u* to 
z* can be accomplished by a single transformation equation, computer run- 


ning time may be as great, or greater, than the table look-up routine, 


1 


due to the number of operations required to calculate the sin ^ or the 


sus 


co 
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APPENDIX V 


FLOW CHART FOR RESOURCE ALLOCATION USING CRITICALITY 
INDEX 


9 


For each T-C plan, form the following computer "word". 


[Activity wo. 


GO TO 
J TH 
ACTIVITY 



















PLACE WOED 
IN Do LIST 


ORDER Al 
WORDS 
BY ETA 


LOOK AT 
NEXT 
LARGER, ETB 








7 ۶ 7 + ۳ 








REMOVE PREVIOD 
WOR? FOR THIS ACT. 
FROM DO LIST 


T=T-AT 









FOR EACH ASTIVIT 
IN Do LIST. CALCULA rar 





RERUN 3 ۱ 
PERT c= 


CALCULATION 
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